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^ ■ Abstract 

The proper incorporation of spin effects in r lepton decays is often of importance. 
In present work the case of the Zj^ — > t + t~ production mechanism is studied in 
^ ' detail. As an example, the effects due to the spin correlations on the potential 

for the Minimal Supersymmetric Standard Model (MSSM) Higgs boson(s) searches 
in the rr decay channel at the Large Hadron Collider (LHC) are discussed. For 
these processes, the Standard Model Z/7* — > rr-pair production is a dominant 
background. The spin effects in high energy physics reactions, can be implemented 
up to certain approximation, independently of the algorithm and matrix elements 
used by the production program. Information stored on every generated event can 
be sufficient. The algorithm based on such approximation is documented. Question 
of the theoretical uncertaintity is partly discussed. 
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1 Introduction 



In a study of "discovery potential" and data analysis of present high energy experiments 
the problems of precise predictions including, simultaneously, signal signatures of the new 
(or studied) physics, backgrounds, as well as all detector related effects should be analysed. 
It is generally believed that a Monte Carlo simulation of the full chain from the beam 
collision to detector response is the most convenient technique to address such question. In 
general it is indispensable to divide Monte Carlo simulation into separate blocks: physics 
event generation and detector response. Later event generation can be divided further 
into parts, describing for example production and decay of the intermediate states. 

In the present paper we will concentrate on the particular class of the processes in- 
volving polarised r leptons. The two main goals of the present paper are: (i) presentation 
of the algorithm for matching r lepton decay and its production, with some control over 
spin effects; in particular in case of Z/7 — ► t + t~ production mechanism, (ii) discussion 
of physical observables sensitive to the spin correlations in the r pair production. 

Spin correlations in the decay of r leptons not only can help to suppress irreducible 
background to the possible resonant r pair production at LHC, such as the MSSM Higgs 
bosons decays, but also help to determine the spin nature of this resonanse. 

In the papers TAUOLA Monte Carlo package for simulation of r lepton decay was 

described. Recently, in Ref. [Q, technical details convenient for using the code in multi- 
purpose environment were collected, and universal interface for combining the simulation 
of r lepton decay, with different packages for generation of physics event was proposed. 
Scheme of Ref. |4[] relies on the information stored in the HEPEVT common block || only, 
and not on the details specific for the production generator, such as PYTHIA (used in 
our examples). In fact, such an interface can be considered as a separate software project, 
to some degree independent both from the specific problem of r production and its decay. 

Our paper is organized as follows: in the next section we will describe new algorithm 
for extracting elementary 2 — > 2 body reaction for // — > Z/7 — > t + t~ , which is necessary 
for properly introducing spin correlations into generation chain. In Sec. 3 we analyze 
spin content of such an elementary function. Sec. 4 is dedicated to the discussion of their 
consequences for the distributions of physics interest. In Sec. 5 we discuss few observables 
where spin effects can improve separation of the Higgs boson signature, in case of the 14 
TeV pp collisions. Summary closes the paper. In Appendix, we explain the basic scheme 
of the spin treatment used in our code. It completes the program manual given in Ref. [|J . 

2 The r polarisation from the Z/7 — > t + t~ decay 

The exact way of calculating spin state of any final state is with the help of the matrix 
element and the rigorous density matrix treatment. This is however not always possible 
or necessary. Often, like in the case of the production and decay of particles in the ultra- 
relativistic limit a simplified approach can be sufficient. Such an approach was developed 
for KORALZ Monte Carlo program |7j and its limitations were studied with the help of 
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Table 1: Probability for the configurations of the longitudinal polarisation of the pair of r 
leptons from different origins. 



Origin 


Pr + 


Pr- 


Probability 


Neutral Higgs bosons: h°, H°, A 


P T+ = +1 
P T+ = -1 


Pr- = "I 
P T - = +1 


0.5 
0.5 


Charged Higgs boson: H + or H~ 


P T+ = +1 


P T - = +1 


1.0 


Charged vector boson: W + or W~ 


P T+ = -1 


Pr- = "I 


1.0 


Neutral vector boson: Z/^y* 


P T + = +1 
P T+ = -1 


Pr- = +1 
Pr- = "I 


Pz 
1-Pz 


Other 


P T+ = +1 
P T+ = -1 


Pr- = +1 
Pr- = -1 


0.5 
0.5 



matrix element calculations of the order a ||. In the following, we study the question 
whwter the approach can be generalised, and the approximate spin correlation calculated 
from the information stored in the HEPEVT common block filled by "any" r production 
program. 

The approximation consists of reconstructing information of the elementary 2 — > 2 
body process e + e~(qq) — > t + t~, buried inside multi-body production process. Let us 
stress that such a procedure can never be fully controlled, as its functioning depends 
on the way the production program fills the HEPEVT common block. It will be always 
responsibility of the user to check if in the particular case the implemented algorithm is 
applicable. Nonetheless our aim is not to replace the matrix element calculations, but 
rather to provide a method of calculating/estimating spin effects in cases when spin effects 
would not be taken care of, at all. Needless to say such an approach is limited (for the 
spin treatment) to the approximation not better than leading-log, and to the longitudinal 
spin degrees only. 

The principle of calculating kinematic variables is simple. The 4-momenta of the 2^2 
body process have to be found. The 4-momenta of the outcoming r's are used directly. 
Initial state momenta are constructed from the incoming and outcoming momenta of the 
particles (or fields) accompanying production of the Z/7 stateQ. We group them accord- 
ingly to fermion number flow, and ambiguous additional particles are grouped (summed) 
into effective quarks to minimise their virtualities. Such an approach is internally consis- 
tent in the case of emission of photons or gluons within the leading log approximation. 

Longitudinal polarisation of r leptons V T depends on the spin quantum number of the 
r mother It is randomly generated as specified in Table [TJ. 



1 The Z/7 state does not need to be explicitly coded in the HEPEVT common block. Note that if 
available, information from the history part of the event, where the 4-momenta of gluons quarks etc. are 
stored, will be used. 

2 The spin quantisation axes are chosen in the same way as in Ref. M. 
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The probability P z used in the generation, is calculated directly from the squares of 
the matrix elements of the Born-level 2^2 process // — > t~t + : 

P _ K / _ T - T+ (+,+) m 

Z \M\ 2 ff -^ T - T+ (+,+) + \M\ 2 ff -^ T+ (-,-) ' 

where / = e, fi,u,d,c, s,b. It can be also expressed (following conventions of Ref. |§), 
with the help of the vector (and axial) couplings of fermions to the 7 ( and Z) bosons. 
Explicit expression for the differential distributions is used. The resulting formula, given 
below, will be useful for the discussion of the numerical results and tests of the next 
sections. 

, x ^p^(s,cos#; 1) 

P ( s 0) = dcos6 \ ' '_>_ (2) 

^f(s,cos6;l) + tP^(s, cos0;-1) KJ 

a cos 6 \ ' > I a cos \ ' ' / 

^^(s,cos#;p) = (l + cos 2 ^)F (s) + 2cos^F 1 (s) 
acosy 

-p[(l + cos 2 9)F 2 (s) + 2 cos9F 3 (s)}. (3) 

with the four form-factors 



Fo(s) = ^(q 2 f q 2 T + 2Re X (s)q f q T v f v T + \x(s)\ 2 (v 2 f + a 2 f )(v 2 T + a 2 T )), 



ira 2 



Fi(s) = ^-( 2Rex(s)q f q T a f a T + \x(s)\ 2 2v f a f 2v T a T ), (4) 



na 2 



F 2 (s) = — ( 2Rex(s)q f q T v f a T + \x(s)\ (v f + a f )2v T a T ) 



ira 2 



Fs(s) = ^ ( 2Rex(s)q f q T a f v T + \x(s)\ 2 2v f a f (v 2 T + a 2 T )) J 



and 



X(S) s-M 2 +isT z /M z - (5) 

Analytic expression for the coupling constants, and their numerical values are given in 
Table g 

In the first step of our discussion the V T is shown as a function of cos#, for several 
centre of mass energies and initial state flavours. The angle 9 denotes r~ scattering angle 
in the Z/^f* rest-frame. It is calculated with respect to the e~, u or d effective beam. The 
Z mass mz = 91.1882 GeV, was taken from Ref. as well as effective sin 2 9w = 0.23147 
and T z = 2.49 GeV. 

In Fig. [I], the angular dependence of the r polarisation for the e + e~ — > t + t~ process 
at the peak of Z resonanse and for the three different values of the centre-of-mass (cms) 
energies above it (upper plot). The r's are strongly polarised in the forward direction, 
where polarisation approaches the value approximately twice as big as the average one 
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Table 2: The 7, Z couplings to fermions. Lowest order approximation, the numerical 
values are given for the effective sin 2 9w = 0.23147. 



Flavour: / 


(If 


v f 




Leptons: f — e, fi,T 


1 


1 + W0 W = 0Q44 
4 sin aw cos 8w 


< . rt /1 — 0.593 

4 sin uw cos c^w 


Quarks : / = u or c 


2/3 


/r| sin2ei ; = 0.227 

4sinoiv costfw 


« ■ r\ n — 0.593 

4sm&w cos 


Quarks : / = d, s or b 


-1/3 


-1+f sm 2 w _ 
4 sin cos 9» 


j ■ /i /i — 0.593 

4 sin tw cos vw 



(the polarisation of Z due to e — e — Z coupling sums with the one due to r — r — Z 
coupling). The polarisation is very small in the backward regions. For cos# = —1 it is 
equal to zero, independently of the centre of mass energy. The reason is the universality 
of both the Z and 7 couplings to all leptons. As can be observed, the r polarisation 
changes significantly (especially in the forward region) with the centre of mass energy. At 
cms energies above the Z-peak, r polarisation is smaller, due to significant contribution 
from the s-channel 7-exchange - production mechanism not contributing directly to the 
polarisation. 

Let us now turn to the production from quarks. In Fig. [I] we show the r polarisation 
for the uu — > t + t~ and dd — > t + t~ elementary processes and the same cms energies 
as in the previous case. For the cms energy close to the Z-peak, the r's produced in 
the forward as well as in the backward directions are strongly polarised. The polarisation 
seem to be almost zero for cos 8 = 0, but in reality is equal (nearly) to the same value as in 
previous case of production from electrons at cos 9 = 0. Also, the average r polarisation 
is close to the case of production from electrons. It is, as it should be, independent 
from the initial state flavour, the Z was produced from. The initial state couplings of 
the Z, affect the angular dependence of the r polarisation only, and give sizable angular 
asymmetry in polarisation. Obviously, the larger the initial state vector couplings to Z 
the larger angular dependence of the polarisation. The above arguments hold, because 
the contribution from the 7 exchange is small. This is the case, in the region of the Z 
peak only. 

Quite different polarisation pattern can be observed for quarks and cms energies far 
above the Z peak (see Fig. |T]). Contribution from the 7 exchange cannot be neglected. 
The 7 — Z interference complicates the pattern even more. Let us comment briefly on the 
numerical results. In the case of u and d initial state, polarisation is negative and nearly 
constant over the forward hemisphere; also the average polarisation is negative. This 
is because of rather large and positive forward-backward asymmetry, and polarisation 
forward-backward asymmetry. The average polarisation increases above the Z-peak and 
approaches, respectively, -35% and -61% in the case of uu and dd annihilation^. Let us 
point, that in case of the r production from the d quarks in the ultra high-energy limit, 
polarisation of the r leptons in the backward directions approaches zero (independently 

3 This may open the way for measuring the flavour of the quarks leading to r pair production. 
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Figure 1: Tests of the TAUOLA universal interface. The r lepton polarisation as a function 
of cos9. We have used ^/~s = M z , 210, 350 and 700 GeV. 

of the numerical value of sin 2 6w). 



3 Towards spin sensitive observables 

Once we have understood the pattern of the r polarisation, let us turn to the question 
of its measurement. We will follow the reasoning similar to the one of Ref. ]TTJ. The r 
polarisation V T can be measured from the energy distribution of its decay products. To 
simplify the discussion of spin effects, the decays r — > ttu were used only[]. We will compare 

4 In case of the 3-body decays, such as r — > e(/j,)i , D sensitivity is smaller and depends on energy of 
leptons; in case of ir — > pv to exploit sensitivity in full, the reconstruction of 7r° is necessary. 
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the cases of vector Z/j* and scalar neutral Higgs boson, produced in hadron collisions and 
decaying into pair of r leptons. The r rest-frame can not be accessed experimentally^, 
the spin effects can be seen, however, through the effects on ir + n~ energy distributions 
and correlations, observed/defined for the laboratory frame. 



5000 - 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.E 



2E„/Vs 




0.1 0.2 0.6 0.4 O.b 0.6 0.7 0.8 0.S 



2E„/Vs 



Figure 2: Single it energy spectrum in the case of r produced from H (left-hand side) or 
Z (right-hand side), yfs = run or -^i = mz respectively. Energies are calculated in the 
rest frame of Higgs boson (or Z). 

Let us start with the following example, where for the production of the r lepton 
pairs Monte Carlo program PYTHIA was used, and for the decay Monte Carlo program 
TAU0LA, and our interface. It was assured, that the invariant mass of the pair of two 
incoming quarks was ^/s = mz = ran- For the time being we discuss energies defined 
in the r + r~ pair rest-frame. With the help of variables z± = 2E n ± / s/s, the spin effects 
are visualized. In Fig. |2] we observe the slope (as expected) of 7r energy spectrum due 
to t polarisation. The slope of the distribution is simply proportional to the polarisation 
(small dip at z± ~ is due to kinematical effect of the it mass). 

^- ~ 1 + V T 2 {z± - 0.5) (6) 

In the case of the plot on the right-hand side of Fig. fj, i. e. decay of the Z, polarisation was 
about -14.7%. In the case of the plot on the left-hand side the spectrum is flat, as would 
be in the case of scalar neutral Higgs boson (or pure 7) where there is no polarisation. 

If polarisation was V T = — 100%, then the distribution slope would be maximally neg- 
ative and touch zero at z± = 1. For V T = 100%, slope would be reversed and distribution 

5 In some cases the invariant mass of Z/j* or Higgs, can be reconstructed from total bilans of the 
observed transverse energies for all tracks in the event. 
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would touch zero at z± = 0. These are the cases of charged Higgs boson and charged W 
boson decays into rv. Respective distributions are shown in Fig. |[ 




G.4 0.5 0.6 0.7 0.8 0.9 




2E„/Vs 



0.4 0.5 0.6 

2E„/Vs 



Figure 3: Single n energy spectrum in the case of r produced from (left-hand side) 
or W ± (right-hand side), y/s = m H ± or yfs = m w ± respectively. Energies are calculated 
in the rest frame of (or W ± ). Continuous line with spin effects included, dotted line 
with spin effects switched off. 



Let us now turn to the question of the spin correlations. As we can see from Table [I], 
the r pairs are produced with the well defined spin configurations ( +, + or — ,— for 
vector bosons; +, — or — , + for neutral Higgs boson), thus, the spin effects are to become 
visible on the two-dimensional distribution build on Z- and z + variables. Indeed, for the 
vector bosons we expect more events when both ir + and ir~ are on the upper side (case 
I) (z + > 0.5, z_ > 0.5 ) or lower side (case II) (z + < 0.5, z_ < 0.5 ), with respect to the 
mixed ones: z + > 0.5, z_ < 0.5 (case III) and z + < 0.5, z_ > 0.5 (case IV). 

The appropriate asymmetry is: 



.4 



FastSlow 



o~i + o~n - o~in - OlV 
o~i + a n + am + aiv 



= 0.24. 



(7) 



This can be compared to the Higgs boson case, where asymmetry A FastS i ow = —0.25. For 
the Higgs boson case the Apastsiow is insensitive on the choice of cms energy, as there is 
no interference effects at all. The sign difference is due to the the vector boson nature of 
Z, as opposed to scalar nature of Higgs boson. 

In order to better visualize the spin correlation effect we have introduced variable 
z s defined as signed part of the following phase space part: surface in z + , z_ variables 
between lines z + = z_ and z + = z_ + a (the sign of the a should be taken). In Fig. |] 
respective plots are given for the H and Z decays. The dashed lines (which in both Z 
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and H cases are flat) correspond to the case when spin correlations are switched off. As 
we expect in the Z/7 — > t~t + decays, due to spin effects a Fast (Slow) 7r ± is most likely 
associated with a Fast (Slow) ir T . The solid line has maximum at z s = and approaches 
zero for z s = ±0.5. Precisely the opposite is true in the case of H — > t~t + . The maximum 
is now at z s = ±0.5 and minimum for z s = 0. 





-0.5 -0-4 -0.3 -0.2 -0.1 



0.2 0.3 0.4 



-0.5 -0.4 -0.3 -0.2 -0.1 



0.2 0.3 0.4 0.5 



Figure 4: Expected number of events as the function of z s variable. Left-hand side plot 
for H , right-hand side for Z . Continuous line with spin effects included, dotted line with 
spin effects switched off. 

Let us now turn our attention to the quantities which (at least in principle) can be 
measured experimentally. Fig. [| shows the invariant mass distribution for the case of 
pure (±, ±) and (— , — ) configurations of the r polarisation and the mass of the resonance 
equal to the Z mass. The spin correlations enhance fraction of the Fast-Fast and Slow- 
Slow configurations which are localised mostly at the shoulders of the 7r + 7r _ invariant 
mass distributions, the Fast-Slow configurations would be localised in the centre of the 
distributions. 

In Fig. |6]we show 7r + 7r~ invariant mass distribution for the Higgs and Z cases. Contin- 
uous line with spin effects included, dotted line with spin effects switched off. Left-hand 
side plot corresponds to the Higgs boson case, right-hand side to the Z. In the case of 
Higgs boson, the mass distribution is peaked centrally, whereas in the case of Z/Y shoul- 
ders of the distributions are more profound, especially the one at lower invariant mass 
additionally enhanced by the polarisation. 

For the vector case the spin correlations enhance fraction of the Fast-Fast and Slow- 
Slow configurations which are localised mostly at the shoulders of the 7r + n~ invariant mass 
distributions. The relative height of the shoulders is sensitive to the average polarisation. 
For the scalar case, the enhancement in the fraction of the Slow- Fast configurations effects 
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Figure 5: The tt + ii invariant mass distribution for y/s = mz- Continuous line for pure 
(— , — ) dotted line for pure (+, +) spin configurations. 




Figure 6: The 7r + 7r~ invariant mass distribution. Left-hand side plot for H ; right-hand 
side for Z/'y*. Continuous line with spin effects included, dotted line with spin effects 
switched off. In the two cases respectively y/s = irtn = rnz- 

in the enhancement of the events localised in the middle of the n + ir~ invariant mass 
distribution, leading to the slightly narrower shape. 
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If all polarisation effects are switched off (dashed lines) the distributions in the two 
cases are identical. That observable, well defined distribution of invariant mass built from 
the visible decay products of the r's, can be helpful in separating Higgs boson signal from 
Z/Y background. 




Figure 7: The 7r + 7r~ invariant mass distribution for uu — > (left-hand plot) and 

dd — > Z/'-f* (right-hand plot) produced with cms energy of 300 GeV. Continuous line with 
spin effects included, dotted line with spin effects switched off. 

The same distribution have been also studied for the off-peak production of Z/j*, i.e. 
for the larger cms energies. In these cases the average polarisation is large and negative, 
also distinct for the uu and dd annihilations. As illustrated in Fig. 0, the effect on the 
7r + 7r~ invariant mass distribution is noticeable. The shape of the distribution might give 
the insight to the structure functions of the colliding protons, once the invariant mass of 
the Z/j* state can be reconstructed. 

4 Case of the Higgs signatures at LHC. 

In the search for new phenomena in accelerator experiments, important parameter to 
estimate discovery chances is the signal sensitivity, equal to the ratio of number of expected 
events from the new physics divided by the square root of the expected background events. 
Any possible increase of such a ratio can improve chances of discovery (or improve limit 
of the exclusion). Use of more sophisticated cuts can be of a great help in a case when it 
can be combined with the physical properties of signal and/or background distributions. 
Equally important is the possibility for the verification of the nature of observed new 
physics, e.g. the quantum numbers of the new resonances. 
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The r leptons are considered as a very promising signature for the searches of the Higgs 



bosons in the Minimal Supersymmetric Standard Model (MSSM) at LHC collider [p~2| . p~3 



Below, we will briefly discuss possible applications of the discussed spin correlations in 
the cases of the neutral Higgs bosons H and A decay into r + r~ pair and the charged 
Higgs boson if ± decays into tv pair. 

The neutral Higgs bosons H and A decays into t + t~ pair, are enhanced for the large 
values of tan f3 (tan (3 denotes the ratio of the vacuum expectation values of the Higgs 
doublets in the MSSM model), with the branching ratio of about 10% for most of the 
range of the interesting Higgs boson mass values ( 150-1000 GeV ). Accessibility of the 
hadronic decay mode of the r + r~ pair has been studied recently by the CMS Collaboration 
Hl4| . The triggering on such events, signal extraction from the irreducible background 
Z/Y — > tt and the reducible backgrounds QCD jets, tt and W + jets, and reconstruction 
of the resonance peak in the r + r~ mass distributions seems feasible for the Higgs boson 
masses roughly above 300 GeV. In that study the r identification is based on the presence 
of a single hard isolated charged hadron in the jet using tracker information. The two hard 
tracks from r~ and r + in the signal events have an opposite sign while no strong charged 
correlation is expected for the QCD jets or W + jets events. The sensitivity of 5a can 
be reached for the large fraction of the MSSM parameter space after three years of data 
collecting at low luminosity. The expected signal-to-background ratio is very large, being 
of the order of one for 5a sensitivity, with the background dominated by the continuum 
Z/-f* — > tt production. The resolution for the reconstruction of the gaussian m TT peak is 
~ 10% of the mass of the Higgs boson. The above performance seems very promising but 
hopefully can be still improved by exploring the spin correlations and polarisations effects. 
The possible improvements may come from the additional suppression of the background. 

But it is also important that exploring effect of the spin correlation on the invariant 
mass of the hadronic decay products might allow to verify hypothesis of the scalar versus 
vector nature of the observed resonance peak in the reconstructed invariant mass of the 
tt pair. As we have seen in the previous sections, the main difference between production 
mechanisms due to Higgs boson versus Z/j* consists of the correlation in energies of the 
r hadronic decay products. 

Let us concentrate on the case of r decays to ttu -most sensitive to the spin correlations. 
In Fig. | the 7r energy spectrum is shown in the laboratory frame for the H (left-hand 
side) and Zj^f (right-hand side) decays. The selection, roughly consistent with what is 
foreseen in the experimental analysis ||12|| , was applied. The minimal transverse momenta 
of the it's were required to be above 15 GeV and the pseudorapidity |?7| < 2.5. Solid 
line shows results with included spin effects, dashed line with the effects switched off. 
Spin correlations lead to the softer spectrum of it in Z/'-f* decays, which may slightly 
suppress the identification efficiency of the r with respect to those produced from the 
Higgs boson decay. The slope in the distribution for the Z/^* (sensitive to the large 
negative average polarisation), depends also on the relative fraction of the uu and dd 
production processes simulated in the proton-proton collision, hence the parametrisation 
of the structure functions. Events were generated with the Monte Carlo PYTHIA 5.7 [13 



and CTEQ2L structure functions; the Higgs boson mass of the 300 GeV and the width 
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Figure 8: The n energy spectrum in the laboratory frame, after basic selection as specified 
in text, for Higgs (left-hand side) and Z/'-f* (right-hand side). Continuous line with spin 
effects included, dotted line with spin effects switched off. 




-G.5 -0.4 -0.3 -0.; 



0.1 0.2 0.3 




0.5 -0.5 -0.4 -0.3 -0.2 -0,1 0.1 0.2 0,3 0.4 0.5 



Figure 9: The expected number of events as a function of the 7r + 7r~ Energy-Energy cor- 
relation variable z s , basic selection as specified in text included. Distributions are shown 
for Higgs (left-hand side) and Z/'-f* (right-hand side). Continuous line with spin effects 
included, dotted line with spin effects switched off. 
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below 1 GeV (as for tan/3 ~ 10) was assumed. For the continuum Z/j*, the cms energy 
of the produced rr pair was taken in the range 300 ±10 GeV. 

Fig. |^ shows the effect of the ttw energy-energy correlations in the cms frame as dis- 
cussed in the previous section. Although the effect seems strong enough to discriminate 
between scalar and vector cases, exploring this effect would require good experimental 
reconstruction of the effective cms frame which might be very difficult. 



The invariant mass distribution of the nir system, see Fig. [10], shows the visible effect of 
including spin-correlations. As discussed in the previous section, in the vector case these 
correlations lead to the more profound shoulders in the distribution, and the relative 
height of them is sensitive to the average polarisation of the produced resonance. In the 
case of the scalar Higgs boson the correlations lead to the narrower distribution, enhancing 
fraction of events localised in its central part. With the expected signal-to-background 
ratio being of one or higher, this effect seems promising for determining spin property of 
the studied resonance. 



N„, 

7000 




Figure 10: The tt + it^ invariant mass distribution after basic selection. On the left for 
Higgs, on the right for Z/^y*. Continuous line with spin effects included, dotted line with 
spin effects switched off. The Higgs mass was assumed to be 300 GeV (see text). 



For completeness let us now turn to the case of the charged Higgs boson. The decay 
into tv pair is a dominant mode below the kinematical threshold of the tb channel, and 
accounts for nearly 100% of cases, almost independently of tan (3. This branching ratio 
is decreasing rather rapidly above tb thereshold, but is still of the order of 10% for the 
Higgs boson mass of 500 GeV and large tan (3. The LHC detectors will be thus sensitive to 



signal ||12l , |14l , with at least 5a significance, for Higgs boson mass up to 400-500 GeV. This 
sensitivity is almost independent on tan (3 below tb threshold. For the Higgs boson mass 
above top-quark mass this sensitivity is expected only for large tan (3. In both cases of 



13 



Higgs boson masses below and above top-quark mass, the main background comes from 
W ± — > r^v decay. Harder pions are expected from the H^ 1 decay than from the W ± 
decays. 

The effect of the spin correlations has been already studied, theoretically in Ref. [IB 
and experimentally in Ref. fll4| . For completeness we show in Fig. [TT] the tt energy 
distribution in the laboratory frame for the decays of H ± and W - spin effects switched 
on and off. The Higgs boson mass of 130 GeV was taken. 




Figure 11: Distribution ofp n /E T after basic selection and in the laboratory frame. On 
the left plot for H ± ( tuh = 130 GeV), on the right for W ± . Continuous line with spin 
effects included, dotted line with spin effects switched off. 
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5 Summary 



We have discussed the spin effects in the r pair production at LHC A few distributions 
presented here, sensitive to the r + r~ spin correlations, can be possibly used for the MSSM 
Higgs boson searches scenarios at LHC to enhance sensitivity of the signal or to verify 
the hypothesis of the spin zero nature of the Higgs boson. 

We have extended the algorithm for the interfacing the r lepton decay package TAUOLA 
with "any" production generator to include effects due to spin in elementary Z/'-f* — > 
t + t~ process. The interface is based exclusively on the information stored in the HEPEVT 
common block. This code is publicly available e.g. from the URL address [ |17|] . 
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Appendix 

Universal interface for TAUOLA package 

Let us recall first some details of the universal interface for TAUOLA and "any" r 
production generator as described in Ref. 0. The interface uses as an input, the HEPEVT 
common block and operates on its content only. As a demonstration example the interface 
is combined with the JETSET generator, however it should work in the same manner with 
the PYTHIA[], HERWIG or ISA JET generators as well. 

The interface acts in the following way: 

• The r lepton should be forced to be stable in the package performing generation of 
the r production. 

• The content of the HEPEVT common block is searched for all r leptons and r neutri- 
nos. 

• It is checked if there are r flavour pairs (two r leptons or r lepton and r neutrino) 
originating from the same mother. 

• The decays of the r flavour pairs are performed with the subroutine TAUOLA. Lon- 
gitudinal spin correlations are generated in the case of the r produced from decay 
of: W —> tv, Z/'-f — > rr, the neutral Higgs boson H — > tt, and the charged Higgs 
boson H — > tv. Parallel or anti-parallel spin configurations are generated, before 
calling on the r decay, and then the decays of 100 % polarised r's are executed. 

6 It was already checked to be the case. 
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In the case of the Higgs boson (for the spin correlations to be generated) the identifier 
of the r mother must be that of the Higgs boson. Here, the particle code convention 
as that used by the PYTHIA 5 . 7 Monte Carlo is adopted. 

In case of the W and Z/7 it is not necessary. If from the same mother as that of the 
r, a v T is also produced, the W is assumed as the mother of the r. Similarly, if from 
the same mother another r with opposite charge is produced, the Z/7 is assumed 
to be the mother of the r pair. 



Photon radiation in the decay is performed with PHOTOS package flq,|19| 



• Let us note that the calculation of the r polarisation created from the Z and/or 
virtual 7 (as function of the direction) represents a rather non-trivial extension. 
Generally, the dedicated study of the production matrix elements of the host gen- 
erator is necessary in every individual case. 

For some of the technical details on how to use the interface we address the reader 
to Ref. 0. For general documentation of TAUOLA to Ref. For the use of PYTHIA 

reference to |15| and references therein will be the best. 



Here, let us concentrate on practical details. The TAUOLA interface is organised in a 
modular form to be used conveniently in "any environments" . 

initialisation 



Initialisation is performed with the CALL TAUOLA (MODE, KEYSPIN), M0DE=-1. All nec- 
essary input is directly coded in subroutine TAUOLA placed in a file tauf ace-jetset . f . 

The following input parameters are set at this call (we omit those which are standard 
input for TAUOLA as defined in its documentation). They are hard-coded in the subroutine. 
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Parameter 


Meaning 


POL 


Internal switch for spin effects in r decy. Normally the user should 
set P0L=1.0, and when P0L=0.0 spin polarisation effects in the 
decays are absent 




(KF=25, 35, 36) Flavour code for h, H and A 


Kr HlOUn 


(KF=37) Flavour code for H + 




yi\r — zoj r lavoui cooe 101 Z/ 


KFGAM 


(KF=22) Flavour code for 7 


KFTAU 


(KF=15) Flavour code for r~ 


KFNUE 


(KF=16) Flavour code for v T 



Event generation 

For every event generated by the production generator, all r leptons will be decayed 
with the single CALL TAUOLA(0,KEYSPIN), (KEYSPIN=l/0 denotes spin effects switched 
on/off). Then, all r leptons, will be first localised, their positions stored in internal 
common block TAUPOS, and the information necessary for calculation of r spin state will 
be read in, from HEPEVT common block. Later spin state for the given r (or r pair) will be 
generated, and finally decay of polarised r will be performed with the standard TAUOLA 
action. In particular, the decay products of r will be boosted to the laboratory frame and 
added to the complete event configuration stored in HEPEVT common block. 

Calculation of the r spin state 

Once CALL TAUOLA (0,KEYSPIN) is executed and r leptons found, the spin states need 
to be calculated. 

First, we look in HEPEVT for the position of r's mothers, and store them, in matrix 
IM0THER(20). Each mother giving r lepton(s) is stored only once, independently of the 
number of produced r's. Later, for every IMOTHER(i) we execute the following steps: 

1. The daughters which are either r leptons, or v T are searched for. 

2. Daughters are combined in pairs, case of more than 1 pair is not expected to be 
important and ad hoc pairing is then performed. 
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3. The two main cases are thus (r r) or (r v). 

4. The default choices are, respectively, Z/7 or W, unless the identifier of IMOTHER(i) 
is explicitely that of neutral (or charged) Higgs boson. 

5. Calculation of the spin parameters is kinematics independent and straightforward 
in all cases except Z/7 (see Section 2 for details on physics). 

6. For Z/7, the P z is calculated with the help of the function PLZAPX (HOPE , IMO , NP1 , NP2) . 
The HOPE is the logical parameter defined in subroutine TAUOLA placed in file tauf ace- 

j etset . f . It tells whether spin effects can be calculated or not. It is set to .false . , 
if available information is incomplete, then, PLZAPX (HOPE, IMO, NP1 ,NP2) returns 
0.5. The IMO denotes position of the r mother in HEPEVT common block NP1 posi- 
tion of r + and NP2 of r~. 

7. To calculate reduced 2 — > 2 body kinematical variables s and cos# subroutine 
ANGULU(PD1,PD2,Q1,Q2,C0STHE) is used. 4-momenta of the incoming effective 
beams and outgoing r + and r~ are denoted by PD1, PD2, Ql, Q2, respectively. 

Run summary 

After the series of events is generated the optional CALL TAUOLA (1 ,KEYSPIN) can be 
executed. The information on the whole sample, such as number of the generated r 
decays, branching ratios calculated from matrix elements etc., will be printed. 

Demonstration program 

Our main program demo.f is stored in subdirectory demo- j etset. It reads in the file 
init . dat which includes some input parameters for the particular run, such as number 
of events to be generated (by JETSET/PYTHIA), or the type of the interaction it should 
use to produce r's, etc. We address the reader directly to the code for more details. It is 
self-explanatory. 
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